(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 7.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     26759,        851]
NotebookOptionsPosition[     26140,        827]
NotebookOutlinePosition[     26478,        842]
CellTagsIndexPosition[     26435,        839]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"p", " ", "=", " ", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"\[Gamma]", "-", "1"}], ")"}], 
    RowBox[{"(", 
     RowBox[{"e", "-", 
      RowBox[{"\[Rho]", " ", 
       RowBox[{
        RowBox[{"u", "^", "2"}], "/", "2"}]}], "-", 
      RowBox[{"\[Rho]", " ", 
       RowBox[{
        RowBox[{"v", "^", "2"}], "/", "2"}]}]}], ")"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"H", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"e", "+", "p"}], ")"}], "/", "\[Rho]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ap", "=", " ", 
   RowBox[{"Sqrt", "[", 
    RowBox[{"\[Gamma]", "  ", 
     RowBox[{"p", "/", "\[Rho]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ea", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"1", ",", "0", ",", "v", ",", "1"}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"u", "-", "ap"}], ",", "0", ",", "uv", ",", 
       RowBox[{"ap", "+", "u"}]}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"v", ",", "1", ",", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"v", "^", "2"}], "-", 
          RowBox[{"u", "^", "2"}]}], ")"}], "/", "2"}], ",", "v"}], "}"}], 
     ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"H", "-", 
        RowBox[{"ap", " ", "u"}]}], ",", "v", ",", "0", ",", 
       RowBox[{"H", "+", 
        RowBox[{"ap", " ", "u"}]}]}], "}"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"eb", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"1", ",", "0", ",", "u", ",", "1"}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"u", ",", "1", ",", 
       FractionBox[
        RowBox[{
         SuperscriptBox["u", "2"], "-", 
         SuperscriptBox["v", "2"]}], 
        RowBox[{"2", " "}]], ",", "u"}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"v", "-", "ap"}], ",", "0", ",", " ", 
       RowBox[{"u", " ", "v"}], ",", 
       RowBox[{"v", "+", "ap"}]}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"H", "-", 
        RowBox[{"ap", " ", "v"}]}], ",", "u", ",", "0", ",", 
       RowBox[{"H", "+", 
        RowBox[{"ap", " ", "v"}]}]}], "}"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"La", " ", "=", " ", 
   RowBox[{"FullSimplify", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"Inverse", "[", "ea", "]"}], "/.", 
      RowBox[{"Solve", "[", 
       RowBox[{
        RowBox[{"ap", "\[Equal]", "a"}], ",", "e"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"\[Rho]", ">", "0"}], ",", 
       RowBox[{"a", ">", "0"}]}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Lb", " ", "=", " ", 
   RowBox[{"FullSimplify", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"Inverse", "[", "eb", "]"}], "/.", 
      RowBox[{"Solve", "[", 
       RowBox[{
        RowBox[{"ap", "\[Equal]", "a"}], ",", "e"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"\[Rho]", ">", "0"}], ",", 
       RowBox[{"a", ">", "0"}]}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"ea", "//", "MatrixForm"}], "\[IndentingNewLine]", 
 RowBox[{"eb", "//", "MatrixForm"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"La", "[", 
   RowBox[{"[", "1", "]"}], "]"}], "//", 
  "MatrixForm"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Lb", "[", 
   RowBox[{"[", "1", "]"}], "]"}], "//", "MatrixForm"}]}], "Input",
 CellChangeTimes->{{3.513760688752114*^9, 3.513760704169449*^9}, {
  3.513761491094468*^9, 3.513761493368145*^9}, {3.513761532497335*^9, 
  3.51376153377159*^9}, {3.513767628370125*^9, 3.513767636241204*^9}, {
  3.513767681107127*^9, 3.513767685202242*^9}, {3.513819022824207*^9, 
  3.513819160513545*^9}, {3.513821973892361*^9, 3.513822063599936*^9}, {
  3.513822161536203*^9, 3.513822217423587*^9}, {3.513822293150993*^9, 
  3.513822386417153*^9}, {3.513822577723125*^9, 3.5138225893488*^9}, {
  3.513830394538936*^9, 3.513830453706557*^9}, {3.513830529877647*^9, 
  3.513830670606396*^9}, {3.513830704329269*^9, 3.513830715460474*^9}, {
  3.513830782641301*^9, 3.513830819921754*^9}, {3.51383202455047*^9, 
  3.513832050903245*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "0", "1", "1"},
     {
      RowBox[{"u", "-", 
       SqrtBox[
        FractionBox[
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]", 
          " ", 
          RowBox[{"(", 
           RowBox[{"e", "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}], 
         "\[Rho]"]]}], "0", "u", 
      RowBox[{"u", "+", 
       SqrtBox[
        FractionBox[
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]", 
          " ", 
          RowBox[{"(", 
           RowBox[{"e", "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}], 
         "\[Rho]"]]}]},
     {"v", "1", 
      FractionBox[
       RowBox[{
        RowBox[{"-", 
         SuperscriptBox["u", "2"]}], "+", 
        SuperscriptBox["v", "2"]}], 
       RowBox[{"2", " ", "v"}]], "v"},
     {
      RowBox[{
       RowBox[{
        RowBox[{"-", "u"}], " ", 
        SqrtBox[
         FractionBox[
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]", 
           " ", 
           RowBox[{"(", 
            RowBox[{"e", "-", 
             FractionBox[
              RowBox[{
               SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
             FractionBox[
              RowBox[{
               SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}], 
          "\[Rho]"]]}], "+", 
       FractionBox[
        RowBox[{"e", "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{"e", "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}]}], 
        "\[Rho]"]}], "v", "0", 
      RowBox[{
       RowBox[{"u", " ", 
        SqrtBox[
         FractionBox[
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]", 
           " ", 
           RowBox[{"(", 
            RowBox[{"e", "-", 
             FractionBox[
              RowBox[{
               SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
             FractionBox[
              RowBox[{
               SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}], 
          "\[Rho]"]]}], "+", 
       FractionBox[
        RowBox[{"e", "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{"e", "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}]}], 
        "\[Rho]"]}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{
  3.513822065836456*^9, 3.51382221837726*^9, 3.513822392238419*^9, 
   3.513822590323016*^9, 3.513830411375941*^9, 3.51383045508492*^9, 
   3.513830587924419*^9, 3.513830631726635*^9, 3.513830673199691*^9, {
   3.513830706988332*^9, 3.513830716616902*^9}, 3.5138308257301*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "0", "1", "1"},
     {"u", "1", 
      FractionBox[
       RowBox[{
        SuperscriptBox["u", "2"], "-", 
        SuperscriptBox["v", "2"]}], 
       RowBox[{"2", " ", "u"}]], "u"},
     {
      RowBox[{"v", "-", 
       SqrtBox[
        FractionBox[
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]", 
          " ", 
          RowBox[{"(", 
           RowBox[{"e", "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}], 
         "\[Rho]"]]}], "0", "v", 
      RowBox[{"v", "+", 
       SqrtBox[
        FractionBox[
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]", 
          " ", 
          RowBox[{"(", 
           RowBox[{"e", "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}], 
         "\[Rho]"]]}]},
     {
      RowBox[{
       RowBox[{
        RowBox[{"-", "v"}], " ", 
        SqrtBox[
         FractionBox[
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]", 
           " ", 
           RowBox[{"(", 
            RowBox[{"e", "-", 
             FractionBox[
              RowBox[{
               SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
             FractionBox[
              RowBox[{
               SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}], 
          "\[Rho]"]]}], "+", 
       FractionBox[
        RowBox[{"e", "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{"e", "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}]}], 
        "\[Rho]"]}], "u", "0", 
      RowBox[{
       RowBox[{"v", " ", 
        SqrtBox[
         FractionBox[
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]", 
           " ", 
           RowBox[{"(", 
            RowBox[{"e", "-", 
             FractionBox[
              RowBox[{
               SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
             FractionBox[
              RowBox[{
               SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}], 
          "\[Rho]"]]}], "+", 
       FractionBox[
        RowBox[{"e", "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{"e", "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
            FractionBox[
             RowBox[{
              SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}]}], 
        "\[Rho]"]}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{
  3.513822065836456*^9, 3.51382221837726*^9, 3.513822392238419*^9, 
   3.513822590323016*^9, 3.513830411375941*^9, 3.51383045508492*^9, 
   3.513830587924419*^9, 3.513830631726635*^9, 3.513830673199691*^9, {
   3.513830706988332*^9, 3.513830716616902*^9}, 3.513830825740513*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", "a", " ", "u"}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], 
       RowBox[{"4", " ", 
        SuperscriptBox["a", "2"]}]], 
      RowBox[{"-", 
       FractionBox[
        RowBox[{"a", "+", 
         RowBox[{"u", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], 
        RowBox[{"2", " ", 
         SuperscriptBox["a", "2"]}]]}], 
      RowBox[{"-", 
       FractionBox[
        RowBox[{"v", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
        RowBox[{"2", " ", 
         SuperscriptBox["a", "2"]}]]}], 
      FractionBox[
       RowBox[{
        RowBox[{"-", "1"}], "+", "\[Gamma]"}], 
       RowBox[{"2", " ", 
        SuperscriptBox["a", "2"]}]]},
     {
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", 
         SuperscriptBox["a", "2"], " ", 
         RowBox[{"(", 
          RowBox[{"u", "-", "v"}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{"u", "+", "v"}], ")"}]}], "-", 
        RowBox[{
         SuperscriptBox[
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], "2"], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], 
       RowBox[{"4", " ", 
        SuperscriptBox["a", "2"], " ", "v"}]], 
      FractionBox[
       RowBox[{"u", " ", 
        RowBox[{"(", 
         RowBox[{
          SuperscriptBox["u", "2"], "+", 
          SuperscriptBox["v", "2"]}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
       RowBox[{"2", " ", 
        SuperscriptBox["a", "2"], " ", "v"}]], 
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", 
         SuperscriptBox["a", "2"]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], 
       RowBox[{"2", " ", 
        SuperscriptBox["a", "2"]}]], 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
        RowBox[{"2", " ", 
         SuperscriptBox["a", "2"], " ", "v"}]]}]},
     {
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", 
         SuperscriptBox["a", "2"]}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], 
       RowBox[{"2", " ", 
        SuperscriptBox["a", "2"]}]], 
      FractionBox[
       RowBox[{"u", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
       SuperscriptBox["a", "2"]], 
      FractionBox[
       RowBox[{"v", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
       SuperscriptBox["a", "2"]], 
      FractionBox[
       RowBox[{"1", "-", "\[Gamma]"}], 
       SuperscriptBox["a", "2"]]},
     {
      FractionBox[
       RowBox[{
        RowBox[{
         RowBox[{"-", "2"}], " ", "a", " ", "u"}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], 
       RowBox[{"4", " ", 
        SuperscriptBox["a", "2"]}]], 
      FractionBox[
       RowBox[{"a", "+", "u", "-", 
        RowBox[{"u", " ", "\[Gamma]"}]}], 
       RowBox[{"2", " ", 
        SuperscriptBox["a", "2"]}]], 
      RowBox[{"-", 
       FractionBox[
        RowBox[{"v", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
        RowBox[{"2", " ", 
         SuperscriptBox["a", "2"]}]]}], 
      FractionBox[
       RowBox[{
        RowBox[{"-", "1"}], "+", "\[Gamma]"}], 
       RowBox[{"2", " ", 
        SuperscriptBox["a", "2"]}]]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{
  3.513822065836456*^9, 3.51382221837726*^9, 3.513822392238419*^9, 
   3.513822590323016*^9, 3.513830411375941*^9, 3.51383045508492*^9, 
   3.513830587924419*^9, 3.513830631726635*^9, 3.513830673199691*^9, {
   3.513830706988332*^9, 3.513830716616902*^9}, 3.513830825836957*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      FractionBox[
       RowBox[{
        RowBox[{"v", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", " ", "a"}], "+", 
           RowBox[{"v", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], ")"}]}], "+", 
        
        RowBox[{
         SuperscriptBox["u", "2"], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], 
       RowBox[{"4", " ", 
        SuperscriptBox["a", "2"]}]], 
      RowBox[{"-", 
       FractionBox[
        RowBox[{"u", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
        RowBox[{"2", " ", 
         SuperscriptBox["a", "2"]}]]}], 
      RowBox[{"-", 
       FractionBox[
        RowBox[{"a", "+", 
         RowBox[{"v", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], 
        RowBox[{"2", " ", 
         SuperscriptBox["a", "2"]}]]}], 
      FractionBox[
       RowBox[{
        RowBox[{"-", "1"}], "+", "\[Gamma]"}], 
       RowBox[{"2", " ", 
        SuperscriptBox["a", "2"]}]]},
     {
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"2", " ", 
          SuperscriptBox["a", "2"], " ", 
          RowBox[{"(", 
           RowBox[{"u", "-", "v"}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{"u", "+", "v"}], ")"}]}], "+", 
         RowBox[{
          SuperscriptBox[
           RowBox[{"(", 
            RowBox[{
             SuperscriptBox["u", "2"], "+", 
             SuperscriptBox["v", "2"]}], ")"}], "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], 
        RowBox[{"4", " ", 
         SuperscriptBox["a", "2"], " ", "u"}]]}], 
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", 
         SuperscriptBox["a", "2"]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], 
       RowBox[{"2", " ", 
        SuperscriptBox["a", "2"]}]], 
      FractionBox[
       RowBox[{"v", " ", 
        RowBox[{"(", 
         RowBox[{
          SuperscriptBox["u", "2"], "+", 
          SuperscriptBox["v", "2"]}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
       RowBox[{"2", " ", 
        SuperscriptBox["a", "2"], " ", "u"}]], 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
        RowBox[{"2", " ", 
         SuperscriptBox["a", "2"], " ", "u"}]]}]},
     {
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", 
         SuperscriptBox["a", "2"]}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], 
       RowBox[{"2", " ", 
        SuperscriptBox["a", "2"]}]], 
      FractionBox[
       RowBox[{"u", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
       SuperscriptBox["a", "2"]], 
      FractionBox[
       RowBox[{"v", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
       SuperscriptBox["a", "2"]], 
      FractionBox[
       RowBox[{"1", "-", "\[Gamma]"}], 
       SuperscriptBox["a", "2"]]},
     {
      FractionBox[
       RowBox[{
        RowBox[{"v", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "2"}], " ", "a"}], "+", 
           RowBox[{"v", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], ")"}]}], "+", 
        
        RowBox[{
         SuperscriptBox["u", "2"], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], 
       RowBox[{"4", " ", 
        SuperscriptBox["a", "2"]}]], 
      RowBox[{"-", 
       FractionBox[
        RowBox[{"u", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
        RowBox[{"2", " ", 
         SuperscriptBox["a", "2"]}]]}], 
      FractionBox[
       RowBox[{"a", "+", "v", "-", 
        RowBox[{"v", " ", "\[Gamma]"}]}], 
       RowBox[{"2", " ", 
        SuperscriptBox["a", "2"]}]], 
      FractionBox[
       RowBox[{
        RowBox[{"-", "1"}], "+", "\[Gamma]"}], 
       RowBox[{"2", " ", 
        SuperscriptBox["a", "2"]}]]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{
  3.513822065836456*^9, 3.51382221837726*^9, 3.513822392238419*^9, 
   3.513822590323016*^9, 3.513830411375941*^9, 3.51383045508492*^9, 
   3.513830587924419*^9, 3.513830631726635*^9, 3.513830673199691*^9, {
   3.513830706988332*^9, 3.513830716616902*^9}, 3.51383082594535*^9}]
}, Open  ]],

Cell[BoxData["H"], "Input",
 CellChangeTimes->{3.513831322951095*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   FractionBox[
    RowBox[{"e", "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
      RowBox[{"(", 
       RowBox[{"e", "-", 
        FractionBox[
         RowBox[{
          SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
        FractionBox[
         RowBox[{
          SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}]}], 
    "\[Rho]"], "/.", 
   RowBox[{"Solve", "[", 
    RowBox[{
     RowBox[{"ap", "\[Equal]", "a"}], ",", "e"}], "]"}]}], "//", 
  "FullSimplify"}]], "Input",
 CellChangeTimes->{{3.513831330982554*^9, 3.513831363717688*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    FractionBox["1", "2"], " ", 
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox["u", "2"], "+", 
      SuperscriptBox["v", "2"]}], ")"}]}], "+", 
   FractionBox[
    SuperscriptBox["a", "2"], 
    RowBox[{
     RowBox[{"-", "1"}], "+", "\[Gamma]"}]]}], "}"}]], "Output",
 CellChangeTimes->{{3.51383135599026*^9, 3.513831364427795*^9}}]
}, Open  ]]
},
WindowSize->{1278, 745},
WindowMargins->{{-7, Automatic}, {-7, Automatic}},
FrontEndVersion->"7.0 for Linux x86 (32-bit) (February 25, 2009)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[567, 22, 4356, 125, 413, "Input"],
Cell[4926, 149, 4306, 131, 183, "Output"],
Cell[9235, 282, 4283, 130, 183, "Output"],
Cell[13521, 414, 5560, 176, 144, "Output"],
Cell[19084, 592, 5867, 187, 142, "Output"]
}, Open  ]],
Cell[24966, 782, 69, 1, 32, "Input"],
Cell[CellGroupData[{
Cell[25060, 787, 668, 22, 60, InheritFromParent],
Cell[25731, 811, 393, 13, 50, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
